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Abstract 

The 1888 paper by Salvatore Pincherle (Professor of Mathematics at the Uni- 
versity of Bologna) on generalized hypergeometric functions is revisited. We 
point out the pioneering contribution of the Italian mathematician towards 
the Mellin-Barnes integrals based on the duality principle between linear dif- 
ferential equations and linear difference equation with rational coefficients. 
By extending the original arguments used by Pincherle, we also show how to 
formally derive the linear differential equation and the Mellin-Barnes integral 
representation of the Meijer G functions. 
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1. Preface 

In Vol. 1, p. 49 of Higher Transcendental Functions of the Bateman Project 
[5j we read "Of all integrals which contain gamma functions in their inte- 
grands the most important ones are the so-called Mellin-Barnes integrals. 
Such integrals were ffist introduced by S. Pincherle, in 1888 [21]; their the- 
ory has been developed in 1910 by H. Mellin (where there are references 
to earlier work) [17] and they were used for a complete integration of the 
hypergeometric differential equation by E.W. Barnes [2]." 



1 



In the classical treatise on Bessel functions by Watson [27|, p. 190, we read 
"By using integrals of a type introduced by Pincherle and Mellin, Barnes has 
obtained representations of Bessel functions 

Salvatore Pincherle (1853 - 1936) was Professor of Mathematics at the Uni- 
versity of Bologna from 1880 to 1928. He retired from the University just 
after the International Congress of Mathematicians that he had organized in 
Bologna, following the invitation received at the previous Congress held in 
Toronto in 1924. He wrote several treatises and lecture notes on Algebra, 
Geometry, Real and Complex Analysis. His main book related to his sci- 
entific activity is entitled "Le Operazioni Distributive e loro Applicazioni 
all'Analisi" ; it was written in collaboration with his assistant. Dr. Ugo 
Amaldi, and was published in 1901 by Zanichelli, Bologna. Pincherle can be 
considered one of the most prominent founders of the Functional Analysis, as 
pointed out by J. Hadamard in his review lecture "Le developpement et le role 
scientifique du Calcul fonctionnel" , given at the Congress of Bologna (1928). 
A description of Pincherle's scientific works requested from him by Mittag- 
Leffler, who was the Editor of Acta Mathematica, appeared (in French) in 
1925 on this prestigious journal [22]. A collection of selected papers (38 
from 247 notes plus 24 treatises) was edited by Unione Matematica Italiana 
(UMI) on the occasion of the centenary of his birth, and published by Cre- 
monese, Roma 1954. Note that S. Pincherle was the first President of UMI, 
from 1922 to 1936. Here we point out that the 1888 paper (in Italian) of S. 
Pincherle on the Generalized Hypergeometric Functions led him to introduce 
the afterwards named Mellin-Barnes integral to represent the solution of a 
generalized hypergeometric differential equation investigated by Goursat in 
1883. Pincherle's priority was explicitly recognized by Mellin and Barnes 
themselves, as reported below. 

In 1907 Barnes, see p. 63 in [T], wrote: "The idea of employing contour inte- 
grals involving gamma functions of the variable in the subject of integration 
appears to be due to Pincherle, whose suggestive paper was the starting point 
of the investigations of Mellin (1895) though the type of contour and its use 
can be traced back to Riemann." In 1910 Mellin, see p. 326fF in [T7j, devoted 
a section (§10: Proof of Theorems of Pincherle) to revisit the original work of 
Pincherle; in particular, he wrote " Before we are going to prove this theorem, 
which is a special case of a more general theorem of Mr. Pincherle, we want 
to describe more closely the lines L over which the integration preferably is 
to be carried out." [free translation from German]. 
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The Mellin-Barnes integrals are the essential tools for treating the two classes 
of higher transcendental functions known as G and H functions, introduced 
by Meijer (1946) and Fox (1961) [6J respectively, so Pincherle can be 
considered their precursor. For an exhaustive treatment of the Mellin-Barnes 
integrals we refer to the recent monograph by Paris and Kaminski [19]. 

The purpose of our paper is to let know the community of scientists inter- 
ested in special functions the pioneering 1888 work by Pincherle, that, in 
the author's intention, was devoted to compare two different generalizations 
of the Gauss hypergeometric function due to Pochhammer and to Goursat. 
Incidentally, for a particular case of the Goursat function, Pincherle used an 
integral representation in the complex plane that in future was adopted by 
Mellin and Barnes for their treatment of the generalized hypergeometric func- 
tions known as pFq{z) . We also intend to show, in the original part of our 
paper, that, by extending the original arguments by Pincherle, we are able 
to provide the Mellin-Barnes integral representation of the transcendental 
functions introduced by Meijer (the so-called G functions). 

The paper is divided as follows. In Section 2 we report the major statements 
and results of the 1888 paper by Pincherle. In Section 3 we show how it is 
possible to originate from these results the Meijer G functions by a proper 
generalization of Pincherle's method. Finally, Section 4 is devoted to the 
conclusions. We find it convenient to reserve an Appendix for recalling some 
basic notions for the generalized hypergeometric functions and the Meijer G 
functions. 

2. The Pochhammer and Goursat generaUzed hyperge- 
ometric functions via Pincherle's arguments 

The 1888 paper by Pincherle is based on what he called "duality principle", 
which relates linear differential equations with rational coefficients to linear 
difference equations with rational coefficients. Let us remind that the sen- 
tence "rational coefficients" means that the coefficients are in general rational 
functions {i.e. ratio between two polynomials) of the independent variable 
and, in particular, polynomials. 

By using this principle with polynomial coefficients, Pincherle showed that 
two generalized hypergeometric functions proposed and investigated respec- 
tively by Pochhammer (1870), see [23], and by Goursat (1883), see [HE], can 
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be obtained and related to each othei0. 

The generalized hypergeometric functions introduced by Pochhammer and 
Goursat considered by Pincherle are solutions of linear differential equations 
of order n with polynomial coefficients, that we report in Appendix. As a 
matter of fact, the duality principle states the correspondence between a lin- 
ear ordinary differential equation {ODE) and a linear finite difference equa- 
tion (FDE). The coefficients of both equations are assumed to be rational 
functions, in particular polynomials. In his analysis [21] Pincherle considered 
the correspondence between the following equations, 

m 

J2 {aho + ahie-' + ah2e-^' + . . . + +ahpe-P') ij^'^^t) = 0, (2.1) 

h=0 

[oofc + aik{x + k) + a2k{x + kf + . . . + amk{x + fc)""] /(x + fc) = 0, 

A:=0 

(2.2) 

where ip{t) and f{x) are analytic functions. These functions are required to 
be related to each other through a Laplace- type transformation il){t) <-> f{x) 
defined by the formulas 

(a) f{x) = m dt , (b) m = ^ e+^* fix) dx , (2.3) 

where / and L are appropriate integration paths in the complex t and x plane, 
respectively. 

The singular points of the ODE are the roots of the polynomial 

ttm + Ctm 1 Z + ttm 2 + • • • + Ctm p = . (2-4) 

whereas the singular points of the FDE are the rots of the polynomial 

aoo + aioz + a2oz'^ + . . . + amoz'^ = . (2.5) 



^In fact, translating from Italian, the author so writes in introducing his paper: "It 
is known that to any linear differential equation with rational coefficients one may let 
correspond a linear difference equation with rational coefficients. In other words, if the 
former equation is given, one can immediately write the latter one and viceversa; further- 
more, from the integral of the one, the integral of the latter can be easily deduced. This 
relationship appears to be originated by a sort of duality principle of which, in this note, 
I want to treat an application concerning generalized hypergeometric functions." 
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For the details of the above correspondence Pincherle refers to the 1885 
fundamental paper by Poincar^ [23], and his own 1886 note [SO]- Here we 
limit ourselves to point out what can be easily seen from a formal comparison 
between the ODE (2.1) and the FDE (2.2). We recognize that the degree 
p of the coefficients in e~* of the ODE provides the order of the FDE, and 
that the order m of the ODE gives the degree in x of the coefficients of the 
FDE. Viceversa, the degree m of the coefficients of the FDE provides the 
order of the ODE, and the order p of the FDE gives the degree in e~* of 
the coefficients of the ODE. 

Pincherle's intention was to apply the above duality principle in order to 
compare the generalized hypergeometric function introduced by Pochham- 
mer and governed by (A. 7) with that by Goursat governed by (A. 6). Using 
his words, he proved that the family of the Pochhammer functions (of arbi- 
trary order p) originates from a linear FDE (of order p) whose coefficients 
are polynomials of the first degree in x , and that the family of the Goursat 
functions (of arbitrary order m) originates from a linear ODE (of order m) 
whose coefficients are polynomials of the first degree m. x = e~* . As a con- 
sequence of the duality principle there is a mutual correspondence between 
the properties of the functions belonging to one family and to the other. 

For the Pochhammer function he started from the ODE of the first order 
(aoo + floie"* + ao2e~^* + . . . + aope"^*) i){t) 

(2.6) 

+ (aio + aiie"* + ai2e"^* + ... + aipe-P*) i)^^\t) = 0, 
to be put in correspondence with the FDE 

(aoo + ciiox) f{x) + [aoi + aii(x + 1)] f{x + 1) + [ao2 + 012(2; + 2)] f{x + 2) 

(2.7) 

+ . . . + [aop + aip(x + p)] f{x + p) = 0. 

In this case Pincherle was able to show that the solution f{x) of the FDE 
(2.7), obtained through the formula (a) in (2.3), depends on p parameters, 
whose logarithms are the singular points of the ODE (2.6). With respect to 
each of these parameters, f{x) satisfies a linear ODE of the Pochhammer 
type of order p . 

^For an account of Poincare's theorem upon which Pincherle based his analysis the 
interested reader can consult the recent book by Elaydi j4j, pp. 320-323. 
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For the Goursat function he started from a FDE of the first order 

^ fix) 

aoi + ail (x + 1) + aai (x + 1)2 + . . . + Omi {x + 1)""] f{x + 1) = 
to be put in correspondence to the hnear ODE of order m 

(aoo + aoie^*) ^(t) + (aio + ane"*) + (020 + 0216"*) ^^^H^) 



Oo0"l"'^102^"l"'j^2 3^^"l"---"l" (^m 

(2.8) 



+ ... + (a^o + a^ie-*) V'('")(t) = 0, 



(2.9) 



Using a result of Melhn, see [HHin], Pincherle wrote the solution of the FDE 
(2.8) as 

where the p^^'s and a^'s are respectively the roots of the algebraic equations 

m 

aoo + aiox + ... + ttmox"^ = a^o Y[{x - p^) = , 

(2.11) 

aoi + ail (x + 1) + . . . + ami (x + 1)™ = ami Ylix - a^) = . 
and c is a constant. If a^O; c^mi are both different from zero, we can assume 

C = — Ctmo/ ^m. 1 ■ 

Pincherle showed that, by setting z = ce*, the ODE of order m (2.9) is 
nothing but the Goursat differential equation (A. 6). 

Furthermore, in the special = , he gave the following relevant 

formula for the solution 

m - ^ p r(.-p )r(.- p.) r(. -,.,.) 

im Ja-ioo i[X — (Ti) i[X — a2) ... 1 (X — (Jm-l) 

where a > ^{pi, p2, ■ ■ ■ , Pm} ■ We recognize in (2.12) the first example in the 
literature of the (afterwards named) Mellin-Barnes integral. 



6 



The convergence of the integral was proved by Pincherle by using his asymp- 
totic formula for T{a + ir]) as 77 — > ±ocE|. So, for a solution of a particular 
case of the Goursat equation, Pincherle provided an integral representation 
that later was adopted by Mellin and Barnes for their treatment of the gen- 
eralized hypergeometric functions pFq{z) . Since then, the merits of Mellin 
and Barnes were so well recognized that their names were attached to the 
integrals of this type; on the other hand, after the 1888 paper (written in 
Italian), Pincherle did not pursue on this topic, so his name was no longer 
related to these integrals and, in this respect, his 1888 paper was practically 
ignored. 

3. The Meijer transcendental function via Pincherle's 
arguments 

In more recent times other families of higher transcendental functions have 
been introduced to generalize the hypergeometric function based on their 
representation by Mellin-Barnes type integrals. We especially refer to the 
so-called G and H functions, briefly recalled in the Appendix. 

In this section (the original part of our paper) we show that by extending the 
original arguments by Pincherle based on the duality principle we are able to 
provide the differential equation and the Mellin-Barnes integral representa- 
tion of the G functions. However, we note that these arguments, being based 
on equations with rational coefficients, do not allow us to treat the Fox H 
functions, since for them an ordinary differential equation cannot be found 
in the general case. 

Our starting point is still the "duality principle" that involves a FDE of 
the first order as in Pincherle's approach for the Goursat function, but, at 
variance of Eq. (2.8), we now allow that the degree of the two polynomial 
coefficients are not necessarily equal. Setting p, q the degrees of these coeffi- 

We also note the priority of Pincherle in obtaining this asymptotic formula, as outlined 
by Mellin, see e.g. [16], pp. 330-331, and [H], p.309. In his 1925 "Notices sur les travaux" 
[22] (p. 56, §16) Pincherle wrote "Une expression asymptotique de r(a;) pour x ^ 00 dans 
le sens imaginaire qui se trouve dans [H] a ete attribuee a d'autres auteurs, mais M. Mellin 
m'en a recemment revendique la priorite." This formula is fundamental to investigate the 
convergence of the Mellin-Barnes integrals, as one can recognize from the detailed analysis 
by Dixon and Ferrar [3], see also [19]. 
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cients, our FDE reads 



aoo + aiox + a2ox'^ + . . . + Upox^ 



+ 



(3.1) 

floi + Oil (x + 1) + (a; + 1)^ + . . . + Oqi (x + 1)^ f{x + 1) = . 

We can prove after some algebra that the associated ODE turns out to be 
independent of the order relation between p and q and reads 

j2 anoi^^'Kt) + e-* ^ a,i V^W(t) = . (3.2) 

h=0 h=0 

As we have learnt from Pincherle's analysis, the solution ilj{t) of the ODE 
(3.2) can be expressed in terms of the solution /(x) of the FDE (3.1), ac- 
cording to the integral representation (b) in Eq. (2.3). 

Now, in view of Mellin's results used by Pincherle (see also Milne-Thomson 
[T8] . §11.2, p. 327), we can write the solution of (3.1) in terms of products 
and fractions of T functions. Denoting by pj (j = 0,1, ... ,p) and cr^ {k = 
0,1, ... ,q) the roots of the algebraic equations 



p 

aoo + aiox + . . . + ttpQxP = Opo LK^ ~ Pj) = ^ ' 

i=i 

<? 

ooi + ail (a; + 1) + . . . + agi {x + 1)*^ = a^i Jlix - (Jk) = . 

k=l 

we can write the required solution as 

111=1 - (Tk) «<?! 

We note, by using the known properties of the Gamma function, that Eq. 
(3.4) can be re- written in the following alternative form 

/(^) = ^^ nrrnt"'"1 . c=(-l)-^+^^. (3.5) 

11^=1 r(l + -x) agi 

Furthermore, introducing the integers m, n such that 0<m<g, < n < p , 
we can combine the previous formulas (3.4)- (3. 5) and obtain the alternative 
form 

^ ' UU+i r(i + Pj - x) nL..+i r(x - a,) ' ^ ^ 
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(3.3) 



(3.7) 



with 

\m+n-p+l ^PO 

We note that Eq. (3.6) reduces to the Pincherle expression (2.10) by setting 
{n — p — q, m — 0}, and to Eqs (3.4), (3.5) by setting {n — p, m — 0}, 
{n — , m — q}, respectively. By adopting the form (3.6)-(3.7) wc have the 
most general expression for f{x) which in its turn allows us to arrive at the 
most general solution ijj{t) of the corresponding ODE (3.2) in the form 



j_ f . nu - pj) nr=i r(i + a,-x 

'^^ ^ 27ri Jl r(i + p,. - x) r{x - 



dx - 



(3.8) 



where L is an appropriate integration path in the complex x plane. 

Now, starting from (3.2) and (3.7)-(3.8) it is not difficult to arrive at the 
general G function namely at its ODE and at its Mellin-Barnes integral rep- 
resentation, both given in Appendix. For this purpose we need only to carry 
out some algebraic manipulations and obvious transformations of variables. 

We first note that using (3.3) the ODE (3.2) reads 



P / d \ _ f d ' 



V^(t) = o. 



Then, putting 

z^ce\ u{z) = ip[t{z)] , ttj fejfc = 1 + (7jfc , 

and using (3.7), we get from (3.9) 



(3.9) 



(3.10) 



:-i)^— .n(.^-a, + i^ 



k=l 



dz 



u{z) = 0, (3.11) 



which is just the ODE satisfied by the Meijer G function of orders m, n,p, g, 
sec (A. 10). Of course, at least formally, the Mellin-Barnes integral represen- 
tation of the G function (A. 8)- (A. 9) is recovered as well and reads (setting 

s = x) 



nr=i m - s) uu r(i - % + 



27ri Jl nLm+i r(i -bk + s) uU+i r(% - 



z" ds. 



(3.12) 
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4. Conclusions 



We have revisited the 1888 paper (in Itahan) by Pincherle on generahzed 
hypergeometric functions, based on the duahty principle between hnear dif- 
ferential equations and linear difference equation with rational coefficients. 
We have pointed out the pioneering contribution of the Italian mathematician 
towards the afterwards named Mellin-Barnes integral representation that he 
was able to provide for a special case of a generalized hypergeometric func- 
tion introduced by Goursat in 1883. By extending his original arguments 
we have shown how to formally derive the ordinary differential equation and 
the Mellin-Barnes integral representation of the G functions introduced by 
Meijer in 1936-1946. So, in principle, Pincherle could have introduced the 
G functions much before Meijer if he had intended to pursue his original 
arguments in this direction. Finally, we like to point out that the so-called 
Mellin-Barnes integrals are an efficient tool to deal with the higher transcen- 
dental functions. In fact, for a pure mathematics view point they facilitate 
the representation of these functions (as formerly indicated by Pincherle), 
and for an applied mathematics view point they can be successfully adopted 
to compute the same functions. In this respect we refer to the recent pa- 
per by Mainardi, Luchko and Pagnini [10], who have computed the solutions 
of diffusion-wave equations of fractional order by using their Mellin-Barnes 
integral representation. 
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Appendix: Some generalizations of the hypergeometric 
functions 

The purpose of this Appendix is to provide a survey of some higher tran- 
scendental functions that have been proposed for generalizing the hypergeo- 
metric function. In particular we shall consider the functions investigated by 
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Pochhammer (1870) and Goursat (1883), that have interested Pincherle in 
his 1888 paper, and the G functions introduced by Meijer (1936-1946), since 
they are re-derived in our present analysis by extending the arguments by 
Pincherle. Our survey is essentially based on the classical handbook of the 
Bateman Project [5] and on the more recent treatise by Kiryakova [9]. 

Let us start by recalling the classical hypergeometric equation. If a homo- 
geneous linear differential equation of the second order has at most three 
singular points we may assume that these are 0, 1, oo . If all these singular 
points are "regular", then the equation can be reduced to the form 

d'^u du 
z{l- z) — + [c-{a + h+ l)z\ — - ahu{z) = , {A.l) 

where a, b, c are arbitrary complex constants. This is the hypergeometric 
equation. Taking c ^ 0,-1, —2, . . . , and defining the Pochhammer symbol 

Y(a + n) 

{a)n = . , i.e. {a)o = 1, , (a)„ = a{a+l) . . . (a+n-l) , = 1, 2, . . . 

T{a) 

then the solution of Eq. (A.l), regular at z = , known as Gauss hypergeo- 
metric function, turns out to be 

u{z) = f:^-f^z-:=F{a,b-c;z). {A.2) 

n=0 ['^)nn. 



The above hypergeometric series can be generalized by introducing p param- 
eters ai, . . .Qp (the numerator-parameters) and q parameters bi, . . . ,bg (the 
denominator-parameters). The ensuing series 

E I?!" 37 := pFqi^i^ . . . , ap; 6i, . . . , 6,; z) , {A.3) 

or, in a more compact form, 

n=0 ^h=l{Oj)n n\ 

is known as the generalized hypergeometric series. In general (excepting 
certain integer values of the parameters for which the series fails to make sense 
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or terminateqj) the series pFg converges for all finite z if p < q , converges 
for \z\ < 1 if p = q + 1 , and diverges for all 2; 7^ if p > g + 1 . The resulting 
generalized hypergeometric function u{z) = pFq will satisfy a generalized 
hyper geometric equation. If we note that Eq. (A.l) satisfied by u{z) = 
F{a,b;c;z) = 2Fi{a,b; c; z) can be written in the equivalent form (see e.g. 
Rainville [25], §46, p. 75) : 




'^+61 ^ ( ^ + 
dz I dz \ dz 



u{z) = 0. 



(A.l') 



we arrive at the equation of order n = g + 1 for u{z) = pFg [(cij)i; (&fc)i; z)] : 



d \ ^ TJ 

j=i \dz ^dz l}^ 





u{z) = 0, 



(A.A) 



The above equation containing the operator zd/dz can be written in a more 
explicit form by using D = d/dz, see e.g. [5] §42, p. 184. Distinguishing 
between the cases p < q and p = g + 1 , we get the following general equations 
in V = v{z) : 



z'^D'^ 
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E 

u=l 



{A^z - B^) D^v + Aov = , p<q, 



{A.5) 



z\\ - + ^ z'^-^ 



{A,z - B,) D-'v + AqV = , p = q + l, {A.6) 



where Aq, A^, B^, are constants. Eq. (A.5) has two singular points, z = 0, 00 
of which 2; = is of regular type, whereas Eq. (A. 6) has three singular points, 
2; = 0, 1, cxD of regular type, hke Eq. (A.l). An equation of the same type as 
Eq. (A. 6) was formerly introduced in 1883 by Goursat [TJ [8] in his essay on 
hypergeometric functions of higher order. 

Another generalization of the Gauss hypergeometric equation was previously 
proposed in 1870 by Pochhammer [23]. He investigated the most general 
homogeneous linear differential equation of the order n {n > 2) of Fuchsian 



"^If at least one of the denominator parameters bk {k — 1, . . . , g) is zero or a negative 
integer, Eq. (A. 3) has no meaning at aU, since the denominator of the general term 
vanishes for a sufficiently large index. If some of the numerator parameters are zero or 
negative integers, then the series terminates and turns into a hypergeometric polynomial. 
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type, namely with only "regular" singular points in {oi, 02, ... , a„, 00} . The 
Pochhammer function thus satisfies a differential equation of the type 



4>ow{z) 







{A.7) 



where the coefficients (puiz) {v = 0,1..., n) are polynomials of degree z/ , 
with (pniz) = (z - ai){z - 02) ... (2: - a„) . 

The pFq functions satisfying Eqs (A.5)-(A.6) and the Pochhammer func- 
tions satisfying Eq. (A.7) are not the only generalizations of the Gauss 
hypergeometric function (A. 2). In 1936 Meijer [12] introduced a new class of 
transcendental functions, the so called G functions, which provide an inter- 
pretation of the symbol pFg when p > q + 1 . Originally, the G function was 
defined in a manner resembling (A. 2). Later [13j, this definition was replaced 
by one in terms of Mellin-Barnes type integrals. The latter definition has the 
advantage that it allows a greater freedom in the relative values of p and q. 
Here, following [5], we shall complete Meijer's definition so as to include all 
values of p and q without placing any (non-trivial) restriction on m and n. 
One defines 



r~iTa,n 



ai, , 



(%)? 



1 

27n 



is) z" ds . 



{A.8) 



where L is a suitably chosen path, 2; 7^ , 2;* := exp [s(ln \z\ + iaigz)] with 
a single valued branch of argz, and the integrand is defined as follows 



is) 



nr=i m - s) uu m 



ttj + S) 



nLr.+ir(i-&fc + s) 



T{aj - s) 



{A.9) 



In (A.9) an empty product is interpreted as 1, the integers 171,11, p,q (known 
as orders of the G function) are such that 0<m<g, 0<'ri<j9, and 
the parameters aj and bk are such that no pole of T{bk — s), k = 1, . . . ,m, 



coincides with any pole of r(l 



«7 + S), J 



,11. 



For the details of the 



integration path, which can be of three different types, we refer to [S] (see 
also [H] where an illustration of what these contours can be like is found). 
One can establish that the Meijer G function u{z) satisfies the linear ordinary 
differential equation of generahzed hypergeometric type, see e.g. p] [p. 316, 
Eq. (A.19)], 



-1) 



^ f d 
j=i\ dz 



ttj + 1 




u{z) = 0. (AlO) 
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For more details on the Meijer function and on the singular points of the 
above differential equation we refer to [S]. Here we limit ourselves to show 
how the generalized hypergeometric function pFq can be expressed in terms 
of a Meijer G function and thus in terms of Mellin-Barnes integral. We have 



1 

271^ 



+ioo 



r(ai + s) 



— a 



(1 



0, (1-6, 



k)l 



■T{ap + s)T{-s) 



T{bi + s)---T{b, + s) ' 
j = l,...,p; |arg(l-2;z) 



ds - 



(All) 



(A.U) 
(A13) 



a,- ^0,-1,-2,... 

Here the path of integration is the imaginary axis (in the complex s-plane) 
which can be deformed, if necessary, in order to separate the poles of T{aj+s), 
j = 1, . . . ,p from those of r(— s) . 

Though the G functions are quite general in nature, there still exist examples 
of special functions, like the Mittag-LefHer and the Wright functions, which 
do not form their particular cases. A more general class which includes those 
functions can be achieved by introducing the Fox H functions [6J, whose 
representation in terms of the Mellin-Barnes integral is a straightforward 
generalization of that for the G functions. For this purpose we need to 
add to the sets of the complex parameters aj and the new sets of positive 



numbers aj and Pk with j = 1, . . 
of (A.8) the kernel G^^^is) into 



,p, k 



q, and modify in the integral 



nLi ^ih - I3ks) U]=i r(l - aj + a,s) 
Lm+i r(l -bk + Pks) ]Tj=n+i r(aj - ajs) 



(A14) 



Then the Fox H function turns out to be defined as 



TTm,n 

p,q 



z 



(%,«i)i=l,. 

{bk, Pk)k=l, 



1 



njm,n , 



s) z^ ds . 



(A15) 



We do not pursue furthermore in our survey: we refer the interested reader 
to the treatises on Fox H functions by Mathai and Saxena pTj, Srivastava, 
Gupta and Goyal [26j and references therein. 
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